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Abstract 

It has been shown that the electron system with ODLRO in the 
reduced density matrix p2 can not support a uniform magnetic field, 
i.e., ODLRO in p2 implies the Meissner effect furthermore, the 

magnetic field trapped in the system is quantized [||. This note extends 
above results in two cases. We show that (1) the system with ODLRO 
in p2 can not support a non uniform, cylindrically symmetric magnetic 
field; (2) the system with ODLRO in p2 can not support a magnetic 
field slowly varying in space, and the magnetic flux trapped in it is 
quantized. 
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Introduction 



In 1962 C. N. Yang proposed that the superconductivity is characterized 
by the existence of off-diagonal long-range order (ODLRO) in the two-particle 
reduced density matrix p 2 for a fermion system 0. However, there was no 
clear-cut proof to show directly that ODLRO in p 2 leads to superconductivity 
for a long time. In recent works it has been proved that ODLRO in 

p2 implies that a non vanishing uniform magnetic field B can not exist in an 
electron system, i.e., there is the Meissner effect. Furthermore, it has been 
shown that the magnetic flux trapped in the system is quantized^. However, 
the question of whether a non uniform magnetic field can be expelled from 
the system with ODLRO in p 2 is still open. The present work extends @ 
and []3| in two cases. We show that 1. the system with ODLRO in p 2 can 
not support a non uniform, cylindrically symmetric magnetic field; 2. the 
system with ODLRO in p 2 can not support a non uniform magnetic field, 
which is slowly varying in space, and the magnetic flux trapped in the system 
is quantized. 

System with Cylindric Symmetry 

Consider the case that N electrons are in a circular cylinder. We chose 
a cylindric coordinates with z-axis along the central line of the cylinder. In 
which r = (R,9,z), and R, 9 and z are unit vectors along the R, 9 and 
z direction, respectively. Suppose that the magnetic field is cylindrically 
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symmetric, and the gauge to be [|] 

A(r) = -6RB(R)R, (1) 

where B(R) is a continuous function. From (1), we have B(r) = V x A(r) = 
B(R)z. Suppose that the system undergoes an infinitesimal rotation 7,89, 

A(r) -> -(6 - 56)RB(R)K = A(r) + V X se(r), 

r (2) 
X se(r) = 56 J R'B(R')dR', 

Ro 

where R is a constant. Eq.(2) shows that the rotation induces a gauge 
transformation. We shall see that it enable us to show that B(r) = as the 
system has ODLRO in the two particle reduced density matrix p 2 . 
The Schrodinger equation of the electrons is 

h ^[f V, + |A(r .)] 2 Vn(ri, r N ) + U(r u r w )^ n (n, r N ) 

(3) 

= £7 n ^ n (n, ...,r N ). 

We assume that the potential U(ri, ...,r N ) depends only on the distances 
between electrons. Under the rotation, (3) becomes 

^EtfV. + fA^.)] 2 ^!,...,^) 

(4) 

+U(t u ...,rjy)^(rx, ...,r N ) = E n ip' n (r u ...,r N ), 

where 

^(r 1; ...,r N ) = il) n (r x - 5r u ...,r N - 5r N ) exp[^| £ X«( r j)]i 

J (5) 

Both { tp n } and are complete sets of orthonormal eigenfunctions of the 
Hamiltonian, which are single- valued in r^. 



The element of the reduced density matrix for two particles can be ex- 
pressed as 

P2(r / i,r' 2 ;r 1 ,r 2 ) =/•••/ ^ - tyl \ ^ exp(-E n /kT) 

x Vn(r'i, r' 2 , r 3 , • • • , r iV )^(r 1 , r 2 , r 3 , • • • , r N ), (6) 

where Z is the partition function, and T the temperature. Substituting ip' n 
for ^ n makes no change for p 2 , since p 2 is basis independent. So we also have 

P2(ri, r 2 ; n, r 2 ) = | • • • J ^—^1 £ exp(-E n /kT) 

x ^n(r'i, r' 2 , r 3 , • • • , r JV )^*(r 1 , r 2 , r 3 , • • • , rjv ), (7) 
Expressing the primed ip' n in (7) in terms of the non primed ip n , we obtain 

p 2 (r'i, r' 2 ; n, r 2 ) = exp{ — [xae(r'i) + X«5e(r' 2 ) - X5<?( r i) - X5e(r 2 )]} 

x^ n (r'i - 5r'i, r' 2 - 5r' 2 , r 3 - 5r 3 , • • • , r N - 5r N ) 
x^*(ri - <Jri, r 2 - 5r 2 , r 3 - 5r 3 , ■■-,r N - Sr N ). 
Shift the integration variables, Q; L — > 0j + 50, « = 3, iV. Note that r; — 5r; = 
(i?,, 0j — 50, Zj). We obtain a relation for p 2 at different space points: 

p 2 (r , 1 , r' 2 ; r 1; r 2 ) = exp{^[xse(r\) + Xse(r' 2 ) - Xsei^i) - X<5e(r 2 )]} ^ 
xp 2 (r' 1 - (Jr'i, r' 2 - 5r' 2 ; r x - 5r x , r 2 - <5r 2 ). 

By definition, the ODLRO in p 2 means that 

p 2 (r' 1 ,r' 2 ;r 1 ,r 2 ) ^ $(r\, r' 2 )$*(r!, r 2 ), (10) 



as |r'j — Tj\,i, j = 1,2, approaches macroscopic distances, and $(ri,r 2 ) 
does not vanish as \r 1 — r 2 | and \r\ — r' 2 | remain finite in microscopic scale. 
Where $ is the eigenfunction of p 2 with the largest eigenvalue. If we assume 
ODLRO in p 2 , (8) and (9) lead to 

$(ri, r 2 ) = exp[ia(£0)] exp{f [xw(ri) + Xse(r 2 )]} 
x$(n -<5ri,r 2 - 5r 2 ), 
where ct(56) is a real function of 59, and it is independent of r 1; r 2 . It is easily 
seen that a(59 + 50') = a(59) + a(<50'), «(°) = °- 

Rotating the system successively, as J2 $@ = 27r, we have J2 = J2 ^ r 2 — 
0. $(ri, r 2 ) is single- valued because ip n is single- valued. Thus, we obtain from 
(10) 

Ri R2 
exp[ia(2n)]exp{^-2n[J RB(R)dR + j RB(R)dR}} = 1. (12) 

-Ro Ro 

(12) leads to 

— a(2vr) + 2vr[ / RB(R)dR + / RB(R)dR] = n— , n = integers. (13) 
Ro Ro 

On the left hand side of (13), R\ and i? 2 can be changed continuously. 
Let R 1 — > i?i + the left hand of (13) gets an infinitesimal increment 
2nR 1 B(R 1 )dR 1 . On the other hand, the right hand side only takes discrete 
values. To be consistent, we have to take B(Ri) = 0. Since Ri is arbitrary 
in the system, we obtain 

B(r) = 0. (14) 

Thus a non vanishing cylindrically symmetric magnetic field can not exist in 
the system with ODLRO in p 2 . In other words, the magnetic field is expelled 
from the system, i.e., there is the Meissner effect. 
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Electrons in the Magnetic Field Slowly Varying in Space 

Consider the case of electrons in a magnetic field which is slowly varying 
in space. The vector potential can be written as 

A(r) = A (r) + Vv?(r), (15) 

and B(r) = V x A (r). In general, <p(r) is multi-valued, however, Vy>(r) is 
single- valued. 

Under an infinitesimal translation Sr, 



A(r) -> A (r - Sr) + Vy?(r - Sr) = A(r) - Sr • V[A (r) + Vy?(r)]. (16) 
In (16) 

Sr • VA = V(Sr -A )-5rxB 

(17) 

= V(Sr ■ A ) - V[Sr • (B x r)] + [(Sr • V)B] x r + 5r x [(r • V)B]. 
Provided that B(r) is slowly varying, such that 

|<?r x B| » \[(5r ■ V)B] x r + 5r x [(r • V)B]| (18) 

i.e., roughly speaking, |VSj| / |B| , i — x,y,z, is much less than L _1 , where 
L is the linear dimension of the system. Under this condition, we can write 

5r ■ VA = V(5r • A ) - V[5r ■ (B x r)], (19) 

and (16) becomes 

A(r) -> A(r) + V[5r ■ (B x r)-5r ■ A - 5r ■ V<p(r)]. (20) 
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We see that the infinitesimal translation induces a gauge transformation. 
In the Hamiltonian of the system, 

H = £ il-AV, + ^A(r,)] 2 + r7(n, r^), (21) 

i 

we assume that the potential U(rx, r N ) depends only on the distances 
between electrons. Suppose that {ip n } is a complete set of single- valued 
orthonormal eigenfunctions of the Hamiltonian. Then {ip' n } is also a complete 
set of single- valued orthonormal eigenfunctions of the Hamiltonian |§, where 

<(r x , r N ) = XsAr > ty^n - 5r, r N - <5r), (22) 

and 

XfrCrj) = Sr-[B( rj ) x iv - A (rj) - V^fo)]. (23) 

Repeating the discussion similar to that leads to (8) ( see also [[§), we 
have 

p 2 (r / i,r' 2 ;r 1 ,r 2 ) = exp{§ [Xfr(r'i) + X<5r(r' 2 ) - Xfr(ri) -X<5r(r 2 )]} 

(24) 

xp 2 (r' 1 - 5r, r 2 - <5r; r x - 5r, r 2 - 5r). 
If there is ODLRO in p 2 , (10) is valid, and 

$(ri, r 2 ) = exp[za(5r)] exp{ — [xst(ti) + X<5r(r 2 )]}$(ri - £r, r 2 - 5r), (25) 

m 

where a(6r) is real, and not dependent on r 1; r 2 . It is easily seen that a(5ri + 
<5r 2 ) = a(5ri) + a(5r 2 ), and a(0) = 0. Repeating successively for infinitesi- 
mal displacements along a closed path, and noting the single-valuedness of 
$(ri,r 2 ), we obtain from (25) the relation 

2ie r 

exp{i^a(<5r)}exp{ — j> dr ■ [B(r) x r - A (r) - Vy?(r)]} = 1, (26) 



where the summation cx(8t) is over the closed path. In a simply connected 

region, around a sufficiently small contour, on which B can be considered as 

a constant vector, we have § dr • [B(r) x r] = 20, § dr ■ A (r) = 0, where 

is the magnetic flux through the surface spanned by the small contour, and 

§ dr ■ V^(r) = 0, £ a(5r) = 0. We obtain, 

2,ie f 2,ie 
exp{— f dr ■ [B(r) x r - A (r)]} = exp{— 0} = 1. (27) 

Thus 

. ch . 

= n -^~i n — integers. (28) 

The contour can be deformed continuously. If B ^ 0, the flux can be 
changed continuously. Which is in contradiction to the right hand side of 
(28). To be consistent, we conclude that B =0. It means that the system 
with ODLRO in pi can not support a magnetic field slowly varying in space, 
i.e., there is the Meissner effect. 

Since B = 0, we can take A = 0. When the path is not simply connected, 
e.g., it winds around an inaccessible tunnel, the path can not be deformed 
across this inaccessible region. In this case we still have J2 cx(8r) = 00, from 
(26), we obtain exp{^ § dr-Vy?(r)]} = 1, thus 

/ch 
dr-V(f(r) = n -^~i n — integers. (29) 

This equation shows that the magnetic flux trapped in the tunnel is quan- 
tized. 

In summary, we have shown that if there is ODLRO in the reduced density 
matrix p2, 1. the system can not support a cylindrically symmetric magnetic 
field; 2. the system can not support a slowly varying magnetic field. In these 
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two cases, we have derived the Meissner effect from the existence of ODLRO 
in p 2 . Furthermore, in the latter case, we have shown that the magnetic flux 
trapped in the system is quantized. We know that the magnetic field is not 
vanishing, e.g., in the boundary region in superconductor. This fact is not in 
contradiction to our result, since the non vanishing magnetic field is in the 
region where it varies rapidly. 
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